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Abstract 

We study the asymptotic behavior of complex discrete evolution equations of Ginzburg- Landau type. 
Depending on the nonlinearity and the data of the problem, we find different dynamical behavior ranging 
from global existence of solutions and global attractors, to blow up in finite time. We provide estimates for 
the blow up time, depending not only on the initial data but also on the size of the lattice. Some of the 
theoretical results, are tested by numerical simulations. 


1 Introduction 

Several theoretical and experimental studies performed on spatially discrete systems, have proved that such 
systems display very reach dynamical behavior, even in the one-dimensional space. The Discrete Complex 
Ginzburg-Landau Equation (DCGL), 

u n - (A + ia){u n -1 - 2u n + u n+ 1 ) = {k + i/3)\u n \ 2 u n + "fUn, (1) 

(where the n-index ranges over the lD-lattice), is a particular discrete evolution equation, whose dynamics can 
lead to extraodinary complicated behavior, ranging from spatiotemporal intermittency and dispersive chaos, 
to self-localization phenomena and the existence of discrete solitons. One of the most interesting applications 
where the DCGL equation 0 may arise, is in the description of the evolution of Taylor and frustrated vortices, 
in hydrodynamic systems of low dimensionality, and it has been proved a fairly reasonable mathematical model 
for investigating phenomena related to weak turbulence EH ESI The famous Discete Nonlinear Schrodinger 
Equation (DNLS) (obtained from 0, in the case A = k = 0) IQ [22] , is encountered in several diverse 
branches of physics, ranging from supeconductivity and nonlinear optics, to the Bose-Einstein condensates. 

The aim of this work, is to provide some results, concerning the asymptotic behavior of solutions of discrete 
evolution equations of the form, 

u n - (A + ia)(u n -i - 2u n + u n+ i) = F(u n )+'yu n , (2) 

^n(0) — U n ^ o, (3) 

where A > 0 and a, 7 £ R, and the classical cubic nonlinearity of 0 , has been replaced by more general 
nonlinear interactions. The lattice may be infinite ( n £ Z) or finite (|n| < N), supplemented with Dirichlet 
boundary conditions. Mainly we are interested in nonlinear interactions of the form F(s) = ( k + i(3)\s\ p (non¬ 
gauge interaction), or F(s) = (k + i/3)|s| p- 1 s, for some p > 1 (gauge interaction). The case a = (3 = 0, 
corresponds to the discrete Ginzburg-Landau equation with real coefficients (DRGL). 

Depending on the type of the nonlinearity, the length of the lattice and the “size” of the initial data, the 
dynamical behavior of solutions ranges from the existence of finite time singularities (blow-up in finite time) to 
the existence of global attractors in appropriate and physically justified phase spaces. 

With respect to the issue of global existence and blow-up of solutions, it has been observed numerically 
that discreteness may have important effects. For example, in the case of the conservative DNLS with gauge 
nonlinear interaction, solutions exist globally, independently of the choice of the initial data and the strength 
of the nonlinearity. This is in contrast with the NLS continuous counterpart, for which solutions may blow-up 
in finite time. As it is shown in Section 3, the DCGL and DRGL equations, serve as a discrete models, whose 
behavior differs with respect to this issue, since both in the case of the gauge and non-gauge interactions, 
solutions may blow-up in finite time, at least in the case of a finite lattice. On the other hand, this behavior is 
in agreement with that of the CGL and RGL partial differential equation (see the recent results of I2ZJ). 


In the case of a dissipative gauge nonlinearity the dynamics of the finite system can be described by a global 
attractor. A simple but interesting observation, is that the entry time to the absorbing ball is independent of 
the size of the initial data, another diferrence with the CGL partial differential equation. To test the theoretical 
estimates for the blow-up time, we perform a first attempt on the derivation of some numerical simulations, for 
the case of non-gauge nonlinearity. 

Motivated by the pioneering work of 0 for lattice dynamical systems of first order and m for extensions 
to various multidimensional lattices of first and second order, we devote Section 4 to the study of in the 

case of an infinite lattice, as an infinite dimensional dynamical system. However, our approach differs from the 
aforementioned references, since we study the DCGL in a particular class of weighted sequence spaces. These 
spaces, cover (but are not limited to) the case of exponential localized solutions for lattice differential equations. 
Substantial physical motivation is given for the study of 0-0 in such spaces, since many of the physically 
interesting solutions of complex discrete evolution equations, present strong spatial localization properties. This 
is true for instance for soliton solutions or breathers 

As in the case of the finite lattice, in the case of a dissipative nonlinearity we prove the existence of a global 
attractor in the weighted phase space. Let us note once again (and in connection with the corresponding result 
for the finite lattice), that discreteness enforces us to impose some restrictions on the parameters, which may 
differ from the usual restrictions on the dissipation parameter, appearing in the case of the continuous CGL 
equation. We conclude by discussing the approximation of the global attractor of spatial localized solutions, by 
the dynamics of the finite dimensional system. 

2 Phase spaces and local existence of solutions 

This preliminary section is devoted to the definition of the appropriate phase spaces and to the basic results on 
local existence of solutions for 0 - 0 . 

The case of the finite lattice for 0-0 (|n| < N), with Dirichlet boundary conditions, will be considered in 
the finite dimensional Hilbert space C 2N+1 endowed with the usual inner product and Euclidean norm, 

n=N n=N 

2 — Re ^ <t>nlfn, IHll : = ^V’GC 2 ^ 1 , (4) 

n=—N n=—N 

We consider now the operators 

{Adlf)\n\<N := '071-1 — 2'(/> n + V’n+l) (B d lf)\n\<N = lfn+1 - ifn- (5) 

It can be easily checked (see also 0 pg. 117]) that 

n=N 

(-Adi/;,ip) 2= ^2 l^n+l - ^nl 2 , (-A d (p,il))2 = (S d 0,S d ^) 2 . (6) 

n=—N 

Hence, we may also consider the inner product and the corresponding norm in C 2Ar+1 , 


n=N 


(0.VOl,2 := {B d (j), B d 1p) 2 + (0,VO2, |M|l,2 := ^ (IV’n+l - tpn \ 2 + IV’nl 2 )- 


(7) 


i=-N 


In our analysis, we shall also use for any 1 < p < oo the norms 


t,=N 


1/p 


= E 


V — I / j I tTL I I 5 imiOO 
n=—N 


= max{|^„| : \n\ < N} , ifeC 2N+1 . 


For any 1 < p < q < oo, there exist constants ci, C 2 depending on N, 

cilHIp < IH| g < c 2 |HI p , g c 2JV+1 . 


( 8 ) 


We note that the norm in 0, is equivalent with the norm in 0 . 


For the case of an infinite lattice (n G Z), a first natural choice for the phase space is to consider complexi- 
fications of the usual real sequence spaces, denoted by 


:= < u = (u„)„ e z <E C : IM \& := ^ |un| p < oo 


V.n(EZ 


Between £ p , spaces the following elementary embedding relation pi pg. 145] holds, 

£ q c£ p , Hk < ||u||* 1 < q < p < oo. 


(9) 


( 10 ) 


For p = 2 we get the usual Hilbert space of square-summable (complex) sequences endowed with the real scalar 
product 


(u,v)p = Re^ u n v n , ii, » £ I 2 . (11) 

n£Z 

Of particular interest is also the existence result in £ x which can be considered as the space of “discrete regu¬ 
larity” , in the sense suggested by H3- 

To cover the situation of spatially localized solutions, we study in weighted spaces, with properly 

chosen weight functions. We consider a weight function 9 n , which is an increasing function of |n|, satisfying for 
all n G Z, the following condition: there exist constants D,d,d> 0, such that 

r i <o n 

(WS) < |0„+1-0„| <D6 n 

[ dOn < 0 n +1 < ddn, 

and we introduce the weighted spaces £^, 

£g = {U n G C : || U || P e p'.= Y^n | U n | P < OO}. 

nGZ 

It can easily be seen that the space £g is a Hilbert space, endowed with the norm || ||^2 and scalar product 

(u, v)e = Re y 0 n u n vu, v G £g. (12) 

nGZ 

Such spaces are the discrete analogue of weighted L p spaces. It follows from (WS), that for any 1 < p < oo, 

£g C £ p , ||u||^p < |H|^, 1 < p < oo. (13) 

Moreover, we observe by using m and (WS), that for any 1 < q < p < oo, similar embedding relations to 
CB9. hold for the weighted sequence spaces, that is 

t-e C £ p e , \\a\\e B < \\u\\e e l<q<p<oo. (14) 

Let us remark that a choice for a weight function satisfying (WS), is the exponential function 9 n = exp(p \ n |) 
for p > 0. Existence of solutions in such spaces will provide us with the existence of (exponentially) localized 
solutions for the DGL equation. An instance where such spaces have been used is in U where the existence 
of exponentially localized solutions has been studied in conservative lattices using a continuation argument, 
related to the anti-integrable limit. 

For local existence os solutions, we shall examine the following examples of nonlinearities F : C —> C: 

(Ni) F(0) = 0 and there exist constants c > 0, p > 1 such that \F(z\) — F(^)| < c(|zi | p_1 + | 22 | p-1 )|al — Z 2 I, 
or alternatively 

(N 2 ) F{z) = /(|z| 2 )^ where /, /' : R — > M, are continuous. 

First, we shall need some information on the nonlinear maps defined by the nonlinear interactions, provided by 
the following 


Lemma 2.1 Let X be either the space C 2N+1 ,£ 2 ,£ g ,£ 1 and assume that F : C —> C satisfies (Ni) or (N 2 ). 
Then the function F defines an operator (still denoted by F) 

F : X —> X , (F(u)) nez := F(u n ), 

which is Lipschitz continuous on bounded sets of X. 

Proof: We focus on the case of the infinite lattice, since the treatment of the finite lattice is almost the same. 
More precisely, we present only the case X = £ g (since the case for £ 2 and £ l is similar). Let u G Br a closed 
ball in £ g , of center 0 and radius R. We have from Q that 

l|F(«)||| < C 2 J2 9 n\Un\ 2P = C 2 \\u\\^ p <C 2 \\u\\ 2 £ ^ (15) 

v ' Q 6 

n£Z 

hence F : £ g —> £ g , is bounded on bounded sets of £ g . 

For u,v G B r, we observe by using the embedding £ g C £°° (provided by (ITol) and O), that 

||F(m) - F(v )\\ 2 2 < <? y; 0 n {\u n \ p ~ 1 + \v n .\ p ~ 1 ) 2 \u n - v n \ 2 (16) 

n(=Z 

< c 2 sup [(|u „| p_1 + K| p_1 ) 2 ] y 0 n \u n -v n \ 2 < c 2 4R 2(p_1) ||m - v\\ 2 2 , (17) 

" eZ raez 6 

justifying that the map F : £ g —> £ g , is Lipschitz continuous on bounded sets of £ g , with Lipschitz constant 
L(R) = c2R(p~ 1 \ 

For the case (N 2 ), we have 


\\ F ( U )\\e e = y 0 n\f(\u n \ 2 )\ 2 \u n \ 2 

n£Z 


( 18 ) 


Since / : K. —* K is continuous, there exists a monotone increasing C 1 -function g : M —> R such that 

\f(p)\ < 5(|p|), for all /)GK, (19) 

(see e.g EH p.g 796]). Now by using flU), we get 

i f(u) h 2 2 = yo n i /(i u n i 2 ) i 2 i u„ i 2 < yd n g{\ u n i 2 ) 2 1 u n i 2 

n£Z n£Z 

< y d n.g (|| U \\%) 2 I u n | 2 < { max g{p)} 2 ye n \ u n | 2 < c(R) || u\\p . 

for some positive constant c(R). Thus we conclude, that the operator F is bounded on bounded sets of £ g . To 
check the Lipschitz property, we may see, that for some v G (0,1), 

\\ F(U) - F(V) \\e s < 2j2dn\f(\u n \ 2 )\ 2 \u n ~V n \ 2 

n£Z 

+2ye n \ f\v I Un | 2 +(1 -v)\v n | 2 ) I 2 (I Un | + | V„ \) 2 \ V n | 2 | U n - V n | 2 • 
n£Z 


For an appropriate C-function gi : M —> K, we get the inequality 


\\F(u)-F(v)\\l 


< 2ye n \ /(|| U || 2 2) | 2 | U n -V n | 2 

n£Z 

+2{ max gi( y p)} 2 c{R) V' 6 » n | it n - u n | 2 

peo, 2 ^ 

< c(R)(\\u-v\\ 2 2g . 


Now in the case of i 1 , we consider as an example the case (Ni). For the Lipschitz condition we have 


||F(u) - F(v)\\ e i < C^(K\ P 1 + \v n \ P 1 )|Wn-'Cn| 

n£Z 

< c sup [(Kl ^ 1 + |u n | p_1 )] ^2 | U n - v n \ < c 2 i?. (p_1) 11 u - v\\(,l. 


n£Z 


n£Z 


This concludes the proof of the Lemma, o 

For local existence of solutions, one could apply alternatively, a semigroup approach or existence theorems of 
ordinary differential equations on Banach spaces, depending on the choice of the phase space. The first approach 
could be of interest, as a starting point for further investigations on the properties of discrete operators m 

2.1 Local existence in £ 2 

In what follows, a complex Hilbert space X 7 endowed with the sesquilinear form Bx{-,-) and the norm 
|| • 11 y , will be considered as a real Hilbert space, endowed with the scalar product (• ■ -)x = ReBx{-,-)- Let 
T : D( T) C X —> X, be a C-linear, self-adjoint, non-positive operator with dense domain D( T), on the Hilbert 
space X, equipped with the scalar product The space Xp, is the completion of D( T) in the norm 

IMIt = I |u||l- — (Tu,u)x for u G X T . We denote by XJ. its dual, and by T*, the extension of T to the dual 
of D (T) , denoted by D(T)*. 

For any u,v € £ 2 we consider the linear operators A , B 7 B* : £ 2 —> £ 2 , 

(Bu') n £.% l = u n _|_i u n7 (.B — u n —i u n7 (20) 

{Au)n£Z — (u n — \ 2 U n T U n -\- 1). (21) 


Lemma 2.2 We assume that £i ,£ 2 > 0. The operator L : £ 2 —> £ 2 , (Lu) ne z = (£1 + i£, 2 )(Au) ne z is the 
generator of a one parameter semigroup of U(t ) on £ 2 , that solves the underlying linear equation it, = L u, 
namely U(t) = exp(Lf). 

Proof: We observe that the operator L, is associated to a non-symmetric bilinear form on £ 2 since for any 
u = u\ + m 2 , v = Vi + iv 2 € £ 2 , 

(Lti, «)p = £1 {(Bui, Bvi)p + (Bu 2 , Bv 2 )p} — £2 {(Bu 2l Bv\)p — {Bu\, Bv 2 )p } . (22) 

However, we observe from that the operator Apu = £1 Au, satisfies the relations 

( A x u,u)i ,2 = —fi\\Bu\\ 2 2 < 0 , (23) 

(Aiu,v)(2 = {u,A\v)(2, (24) 

therefore defines a self-adjoint operator on D(A) = X = i 2 and A\ < 0. We denote next by £ 2 ^, the Hilbert 
space with the following scalar product and induced norm 

(u,v)i 2 := ^ 1 (Bu,Bv)p + (u,v)p, (25) 

M% ■= Hi\\Bu\\l + \\u\\l. (26) 

The usual norm of £ 2 and (|26l) are equivalent (see also J 22 j), since for some constant c(£i) > 0, 

I Ml?* < \[u\\% <c\\u\\ 2 2 . (27) 

Note that the graph norm 

IMI-DLU) = ll^l u ll ? 2 + |M|? a > 

is an equivalent norm with the f-norm since 

I Ml ? 2 < Cl ^2 \ Un + 1 ~ 2 u n + u n - 1| 2 + \u n \ 2 < c||u|| 2 2 . 

n£Z n£Z 


In our case, as it is indicated by 03)-0Il>-<E3, we may choose Xa 1 = f 2 , equipped with the norm HuH^ = 
||m||^ — ( Aiu,u)x = ||m ||^2 , for u £ fi 2 . Moreover, D(Ai ) = X = l 2 = D(Ai)*. Obviously A\ = Ai and A\ is 

the generator of a strongly continuous semigroup on l 2 . 

Thus, we may consider L, as a perturbation of A\ by the bounded (skew-adjoint) linear operator A 2 : i 2 —> fi 2 , 
(A 2 u) ne z = i&(Au) ne z, and apply [2B1 Theorem 1.1], to justify that L = Ax + A 2 , is the generator of a strongly 
continuous emigroup on l 2 . o 

With Lemma IP at hand, and applying it for the case £1 = A, £2 = a, we may recast 0-0 into the form 
of the integral equation 

u(t) = U(t)uo + f U(t — s)Fi(u(s))ds , Fi(u) := — ( F(u ) +7 u). (28) 

Jo 

With the help of Lemma IP we can handle 1281) by a contraction method and the local existence result can be 
stated as follows (we refer to |TTII.|T5] for the proof). 

Theorem 2.1 We assume that A, a > 0, and conditions (Ni) or (N 2 ) are satisfied. Then there exists a function 
T* : l 2 —> (0, 00 ] with the following properties: 

(a) For all uq £ fi 2 , there exists u £ C([0, T*(uo)), fi 2 ), such that for all 0 < T < T*(uo), u is the unique solution 
of 01-01 in C([0,T],£ 2 ) (well posedeness). 

(b) For all t £ [0,T*(tto)), 

T *(u 0 )-t> ^ :=T R , R = 2\\u(t)\\p, (29) 

where L(R), is the Lipschitz constant for the map F\ : fi 2 —> £ 2 . Moreover the following alternative holds: (i) 
T*(u 0 ) = 00 , or (ii) T*(uq ) < 00 and lim t i- 2 '»(„ 0 ) ||u(t )||^2 = 00 (maximality). 

(c) T* : £ 2 —» (0, 00 ] is lower semicontinuous. In addition, if {unojneN is a sequence in £ 2 such that u n 0 —> uq 
and if T < T*(uq), then S(t)uo n —■► S(t)uo in C([0,T],£ 2 ), where S(t)uo = u(t), t £ [0, T*(uo)), denotes the 
solution operator (continuous dependence on initial data). 

We note that Lemma IP and Theorem IP remains valid in the cases A = 0, a > 0 (DNLS type equation) and 
A > 0, a = 0. 

2.2 Local existence in K. 2 e and i 1 

Since the operator A is not symmetric in £g, we cannot apply the analysis for the operator L and Theorem 12 .II 
for the local existence in £g. On the other hand for the case of I 1 , Hilbert space methods are not applicable. In 
both cases, the problem can be treated by general existence Theorems in Banach spaces. 

Concerning the operator L : X —> X, X = fig, I 1 we have the following 

Lemma 2.3 The operator L : A' —> X, defined by is globally Lipschitz on X. 

Proof: Let u, v £ B R . Then (Lu)„ e z - (Lu) rieZ = (A + ia) {(u n+1 - fn+i) - 2(u n - v n ) + (« n _ 1 - u„_i)} and 

it follows that || Au — Av\\x < L\\u — v\\x, where L is independent of R. o 

Thus,the local existence result in the case of the spaces A' = fig,! 1 , can be stated as follows. 

Theorem 2.2 We assume that A, a > 0, and conditions (Ni) or (N 2 ) are satisfied. For all uq £ X, there 

exists T*(u 0 ) > 0, such that for all 0 < T < T*(ufi), there exists a unique solution of the problem 0-0, 
u(t) £ C 1 ([0,T],.X'). 

Proof: This time, we write 0-0, as an ordinary differential equation in fig 

u(t) + $(u(i)) = 0, 

«(0) = u 0 , 

where <h(n) = F(u) + L(u) + *yu, and u(t) lies in A. Lemmas 12.112.31 suffice for the application of standard 
existence and uniqueness Theorems for ordinary differential equations in Banach spaces (e.g. generalized Peano 
and Picard-Lindelof Theorems) (3S, pg. 78-82], o 

We remark that IP holds also in the case where A > 0 and a, 7 £ R. 





3 The case of a finite lattice with Dirichlet boundary conditions 

We will study in this section the asymptotic behavior of solutions to the DCGL and DRGL equations, considered 
in a finite lattice, assuming Dirichlet boundary conditions 


ia)(u n - i 

- 2 U n + U n+ 1) = F(u n ) + 7 U n , 1 77,| < N, 

(30) 

-(JV+1)(’) 

= U (N+1){‘) = 

(31) 

u n (0) 

= «n,o, \n\ < N. 

(32) 


Our questions concerning the life span of solutions to m-m consider nonlinear interactions of the following 
forms 

(A fQ) F(s ) = (k + i(3)\s\ p for some p > 1. (non-gauge type nonlinearity), 

(0) F(s) = (k + */3)|s| p_1 s for some p > 1 (gauge type nonlinearity). 

First we present theoretical estimates for the blow-up time, for several physically interesting parameter 
regimes and we conclude with the existence of a global attractor in the case of the gauge nonlinearity. 


3.1 Blow-up in finite time for the case of a DCGL equation in the case of non¬ 
gauge nonlinearity 

A. (A > 0 and (3 > 0). Motivated by [37], for any t £ (0 ,T*) we define the function 

n=N 

M(t) = Im y «„(*), L = 2N+1, fj > 0, (33) 

n=—N 


and we assume that Im u o,n > 0- The unspecified parameter ex, will be related to a scaling argument, 

on the investigation of the behavior of the upper bound for the blow-up time, that we shall derive in the sequel 
(see Remark EHt . 

We differentiate ifTIill to obtain 


M'(t ) 



n=N 

Y. (-JUn + U n ) 
n=—N 


(3e ~^ 

L a 


n=N 



K| p > o. 


e 7 ‘ 


L a 


n=N 

Im {(A + ia)(u n -i 

n=—N 


— 2 u n + u n + 1 ) + F(u n )} 


(34) 


Now, an application of inequality 0 to m - implies that 


M{t) < 



i/p 


with q = —r, hence 

^ p— 1 ’ 


n=N 

M(t) p < e-P^LP - 1 - 1 ” 7 Y \u n \P. 

n=—N 

Inserting iBHll to iTHl) . we derive the inequality 

M'(t) > 

Now using 1,361) . and differentiating the function M 1 p (t) we observe that 

— (MUy-P) < -(p- l)^e (p - lht L (1 - p)( 1 ~ a) . 
at 


(35) 


(36) 


(37) 








Integration of m with respect to time, implies that 

M( 0 ) 1-p - | ( e (P- 1 b* - 1 ) 7 ^ o, 

M(0) 1_p - (p - 7 = 0 . 


Since M(t) > 0 for all t G [0, T*), we obtain from iTTHll . that the maximal existence time T* can be estimated as 

(38) 


T < 

Note that we have assumed that 


In {1 + lM (O) 1 -^-* 1 -^ 1 -")} , 7 ^ 0 , 

0) 1 -* 7 = 0. 


^M(0) 1_p > 7^0. 

P 

B. (A > 0 and k > 0) This time, we consider the quantity 

—yt n=N 


N(t) = ~— Re ^2 «*(*)> L = 2N + 1, 


(39) 


(40) 


%=-N 


assuming now that R eJ2n=-N w o ,n > 0. We observe that 


%=N 


Up- 7* w _ 

N'(t) = ^— J2 l«n| p >0. 


n=—N 


Following similar arguments to those we used in case A., we obtain that the maximal existence time T* can be 
estimated as 


T* < 


™ In {1 + 2 7V (0) 1 - p L-( 1 - p )( 1 -)} 

l -p Lr} j-pKi-o-) 


iV(O ) 1 


(p-i)fc 


7^0, 

7 = 0. 


(41) 


This time we have assumed that 

^ N ( 0) 1 -* > _x(i-p)(i-«0 7^0. 
k 


(42) 


We summarize the above results, in the following 

Theorem 3.1 A. We assume that X > 0, /? > 0 and Im X)"=-iv u o,n > 0. T/ien for the DCGL equation 
with nonlinear interaction (AfQ), the maximal existence time is estimated by I.Vfll ) . 

B. We assume that X > 0, k > 0 and Re u 0 -ji > 0. TTien /or t/ie DCGL equation f 3(A) - 1321) with 

nonlinear interaction {AfQ), the maximal existence time is estimated by m - 

Remark 3.1 (Scaling limit and Blow-up). The estimates for the upper bound of the blow up time T*, 
can be interpreted in the following way, employing some scaling arguments. The parameter a is unspecified in 
the above argument. Let us consider the behaviour of the upper bound for the blow up time T*, as we take 
the case of a large system (L —> oo). Since p > 1, and under the assumption that M{ 0) = 0(1), we see that 
the upper bound tends to zero if a > 1, whereas the upper bound tends to infinity if a < 1. The upper bound 
is independent of L if a = 1. This means that if a > 1, then in the limit of large systems L —> oo, the system 
blows up instantly, when a = 1 , one may ask if the blow-up time is independent of the lattice size L , whereas 
if a < 1 , one may ask if the system may have longer lifetimes. 

The above observation, along with the condition that M( 0) = 0(1), allows for some heuristic investigations, 
regarding the lifetime of solutions and its possible dependence on the way the initial data decay: from the 
definition of M( 0), we see that this quantity is of order 1, as long as the sum of u„(0) on all lattice sites, scales 
as L a . That means that if u n { 0) ~ n s as n —* oo, then M(0) ~ L 1+5 ~ a , so that M(0) = 0(1), as long as 
S = a — 1. Thus, we may have instant blow up, as long as <5 > 0. On the other hand, it seems to be an 
interesting question, if the solution may live for longer times as long as 5 < 0 (since the behavior of the upper 
bound does not necessarily imply a similar behavior of the blow-up time). Finally, it appears that the upper 
bound is independent of the size of the system, as long as d = 0. In conclusion, if the initial data do not decay 
fast enough in space, the solution will blow up instantly for large systems, whereas for spatially decaying initial 
data, one may conjecture that the solution may live for longer finite times . 1 

1 However, in practice the actual life time may be smaller, as it can be seen from numerical simulations. 









Remark 3.2 (Indications for global existence). We observe that conditions ifTiili and 01 are always valid 
in the case where 7 > 0 (i.e in the case where the linear term acts as a linear source). On the other hand in the 
case 7 < 0 (linear dissipation), these conditions Q 1 and imply that 


/,(!—P)(! —o') 

°<- 7 < . /5 > 0 , 


0 < -7 < 


kL (l-p)(l-a) 


k > 0 , 


iV(0) 1-p 

providing a range for the dissipation parameter, for a possible observation of blow-up in finite time. 


3.2 Blow-up in finite time for the case of a DRGL equation in the case of gauge 
nonlinearity 

We shall examine now the case of a DRGL equation [a = 0 and (3 = 0) in the case of the gauge type nonlinearity 
(Q) and k > 0. For this case, we consider the scalar quantity 


E(u) = 



(BdU)n l 2 -^ E 


n—N 


U n - 


i,=-N 


p + 1 


n=N 

E 1 

n=—N 


|p+l 


and now we shall consider the function 


M(i) 


1 

L a 


n—N 

E MOI 2 - 

n=—N 


(43) 


(44) 


Now we multiply equation m in the C 2Ar+1 scalar product. We get the energy equation 

, n=N n=N n=N n=N 

2 j t E M 2 + A E \( B d u )n\ 2 -7 E M 2 - fc E M p+1 = 0- 

n=—N n=—N n=—N n=—N 


Then from 01 and 01 we obtain 


2A 


i=JV 


n=—N 


«'<(> = E + Y. 


2=N 


\Un 


= -4E(»(i)) E 


n=—N 
n=N 


Ok n = N 

2 + # E i 

n=—N 


IP+1 


|P+1 


(p + 1 )L° 


i=-JV 


(45) 


(46) 


Multiplying the DRGL equation by and keeping real parts, we observe that E{u{t)) < E(uq). Now under 

the assumption that E(uq) < 0, we get from iHhl) that 


M'(t) > -4E(u 0 ) + 


2 k{p — 1 ) 
(p + 1 )L a 


n=N 

E ki p+1 

n=-AT 


> 


2 k(p — 1 ) 
(p+ l)L° 


n=N 


E ki p+1 - 

n=—N 


(47) 


Once again, inequality © implies that 


n=TV / n=N 

E \ U n\ 2 < L ^{ E K ! P+1 

n=—N \n——N 

Therefore, we have for M (t) that 

n=N 

M(tr<L e i«.r, 

n=—N 



(48) 


(49) 










(50) 


Now setting p\ = we insert into lITHl and we obtain 


9 k 

M\t) > 

Pi 

We proceed as for the derivation of the estimate EH- Since p > 1, we have from pH) that 

at pi 

and integration with respect to time, implies that 

9k 

M(t) 1 ~ p < Mfoy-P -(p- 1))—L^-pX 1 -")*. 

Pi 

Since Af(f) > 0 for all t £ [0,T*), we have that this time, the maximal existence time T* satisfies the estimate 


r>< (P + 1) L-V-M 1 -") 
~ k(p- l) 2 M(oM 


(51) 


Summarizing, in the case of the gauge nonlinearity, we have the following 

Theorem 3.2 We assume that A > 0, k > 0 and that the initial energy is such that E(uq) < 0. Then for the 
DRGL equation m m (a = P = 0) with nonlinear interaction (Q), the maximal existence time is estimated 

by EH). 


Remark 3.3 (Scaling limit and Blow-up) Similar comments as in Remark 13.11 hold with respect to scaling of 
the initial data and blow-up, for the case of the DRGL equation I®)-(EH, with the gauge nonlinear interaction 

(S). 


3.3 Comparison of DRGL equations with DNLS, with respect to global existence 
of solutions and blow-up in finite time 

The fact that solutions of DRGL, in the case of gauge nonlinearity (G), may blow-up in finite time, is in 
contrast with the behavior of solutions of the discrete nonlinear Schrodinger equation (DNLS) with the same 
nonlinearity. The solutions of DNLS in the case of gauge nonlinearity, exist globally, unconditionally with 
respect to the degree of the nonlinearity, the size of the initial data and the sign of the initial energy as it 
was observed first numerically in |5J. A detailed discussion on the asymptotic behavior of solutions of DNLS 
equations is presented P2jj- For the sake of completeness and for a comparison, we present here the simple proof 
on global existence of solutions, in the case of the infinite lattice. That is, we shall consider the DNLS lattice 
differential equation (the case A = fc = Oof0-@) 

u n -ia(u n _i - 2u n + u n+1 ) = iP\u n \ p ~ 1 u n , n £ Z, 1 < p < oo, (52) 

it„(0) = u n} o, n £Z. (53) 


and we assume that a, P > 0. Theorem rm (or m covers also the case of DNLS (EH-®): For all u 0 £ £ 2 , 
there exists u £ C([0, T*(u 0 )), £ 2 ) such that for all 0 < T < T*(u 0 ), u is the unique solution of EH-EH in 
C([0,T],^ 2 ). Taking the scalar product of EH) with iu, we obtain that 


^11^(011^=°, or Il u (0llf2 = \\uo\\%, for every t £ [0,T*(wo)), (54) 

Although the conserved quantity EH suffices to demonstrate global existence, to elucidate the interplay of 
nonlinearity and discreteness, we shall examine the DNLS Hamiltonian 

E M0) = EM, E( U (t)) :=f£ \{Bu) n {t)\ 2 -£ K(0l P+1 - (55) 


From EH and EH we may easily derive the ’’conservation law” 

E 1 («(0)=E 1 (« 0 ), Ei(ti(0) := ||N0ll! 


P 


p+ 1 


u (0ll£+i> 


(56) 






Figure 1: Theoretical upper bound for the blow-up time (T*) and numerically obtained (T s i m J blow-up time for 
the finite lattice, in the case of initial data for which a < 1 . 


where the £f is defined by for £1 = 1. Then, by using m- <Gni) and we may derive the estimate 

Mt )|| < IK||| + ^^{|K||^ + IK 0 II^} 

< IKI| + ^^{lKII? 2 +1 + IKi)ll? 2 +1 } 

< IKII + ^^ylKir 1 . (57) 

As a consequence of E3 we obtain that T*(uq) = oo and sup |||u(t)||^ 2 , t £ [0, oo)| < oo. The proof is very 
similar in the case of the finite lattice, assuming Dirichlet boundary conditions. 

Let us mention that this behavior of the DLNS system, is not only in contrast with the DRGL system, but 
also with its continuous counterpart, 

d t u — iau xx = i/3\u\ p ~ 1 u, x £ M, t > 0, (58) 

u(x,0) = uq(x). 

In order to clarify these differences, let us recall the main results concerning (I5%1) (see fTUlimi?5| l: For u o £ H 1 (R) 
and 1 < p < oo there exists a unique maximal solution of iJSSJ), u{t) £ C([0, T max ), H 1 (R))nC 1 ([0, T max ), L 2 (R)). 
In addition: If 1 < p < 5 then T max = oo and u is bounded in H 1 (R). Let p > 5. Assume that uq £ H 1 (R) 
such that f R |a;| 2 |ito| 2 da" < oo (initial data with finite variance) and E(uq ) < oo. Then T max < oo. On the other 
hand if HuoIIh 1 is sufficiently small, T max = oo and u is bounded in H 1 (R). 

In the case of DNLS, the assumption of ininitial data with finite variance, reads as Ylnez l n l 2 l M n,,o | 2 < °°; 
and such data belong to £ 2 . On the other hand, it follows from Theorem im A., that the solutions of DNLS 
system with non-gauge nonlinear interaction ( MQ) may blow up in finite time (at least in the case of Dirichlet 
boundary conditions). Thus, regarding the DNLS system, with respect to global existence of solutions and 
blow-up in finite time we may comment with the following 

Theorem 3.3 A. Consider the DNLS equation I5P with Dirichlet boundary conditions and nonlinear inter¬ 
action (A/’t/). Assume that a, f3 > 0 and Im Yln=-N u o,n > 0. Then the solution blows-up in finite time for all 






























































Figure 2: Theoretical upper bound for the blow-up time (T*), and numerically obtained (T s - im ) blow-up time for 
the finite lattice in the case of initial data for which a = 1 . 


1 < p < oo. The maximal existence time is estimated by IM) . 

B. Consider the DNLS equation IM) with Dirichlet boundary conditions and nonlinear interaction ( Q ). Assume 
that a, [3 > 0. Then the solution exists globally in time unconditionally with respect to the initial data and the 
sign of the initial energy, for all 1 < p < oo. The same holds for the DNLS infinite lattice That is for any 
«o € f 2 , the solution of IM) is in C 1 ([0, oo), i 2 ). 

3.4 Numerical simulations 

Eventhough the aim of the paper, was to show -using some analytical arguments-fmite time blow-up of solutions 
of discrete complex lattices, we decided to test the theoretical estimates for the upper bound of the blow-up 
time, and the heuristic scaling arguments of Remark EH1 numerically, against the observed blow-up times, for 
some parameter values. 

The finite lattice equations, have been integrated numerically over time, using a fourth order Runge Kutta 
scheme, implementing Dirichlet boundary conditions. For the case of non-gauge type nonlinearity, and for 
initial data with cr < 1 , the numerically obtained blow up time T s im, is shown in figurc^and compared with the 
theoretical estimate for the upper bound of T* (still denoted for simplicity, by T*). 

In the first graph, we show the variation of T s j m and the upper bound for T*, as a function of the number of 
lattice sites N, in the second with respect to 7 , in the third with respect to p and in the fourth with respect to 
p. In the first graph, it seems that the observed blow-up time is independent of N (although inspection of the 
data, show a slow increase). Note however, that according to (l.'lHIl . and as clearly stated in Remark Id. II when 
cr < 1, the behavior of the upper bound, does not necessarily imply a similar behavior of the blow-up time. In 
figure |2j we demonstrate the case of initial data with cr = 1 . We first observe, that the scaling limit argument 
we propose, with respect to the number of lattice sites, and the theoretical estimates for the variation of the 
upper bound ,with respect to the various parameters discussed above, seem to capture-at least quatitatively- 
the variation of the numerically observed blow-up times, with respect to these parameters: in the case a = 1 , 
the numerically obtained blow-up time seem to be independent of the lattice size L = 2TV + 1. Furthermore, 
regarding the variation with respect to the parameters (and especially with parameter /?), it is surprising that 
the upper bound is close (same order of magnitude) with the numerically observed blow-up times. We also 
observe that the trend predicted by the theory seem to be verified, as far as the variation with respect to the 
parameters is concerned. 























































Figure 3: Theoretical upper bound for the blow-up time (T*), and numerically obtained (T s - im ) blow-up time for 
the finite lattice in the case of initial data for which a > 1 . 


Finally in figure|3]we repeat the same calculations with initial data with a > 1. Especially in the case of the 
variation with respect to the number of the lattice sites, the numerical simulations seem to verify the instant 
blow-up of solutions for increasing lattice size, as it is predicted by the proposed scaling argument. 


3.5 A global attractor for the case of non-gauge nonlinearity 

We conclude our study for the finite complex lattice, by verifying existence of a global attractor, for the case of 
a dissipative gauge nonlinearity. 

Proposition 3.1 We assume that A > 0, a,/3, 7 G R. and k < 0. Let uq := (u ra ,o)|n|<iv € C 2JV+1 . For UW - 
ViiA) . with nonlinear interaction (G), there exists a unique solution \:i(A ) - l,VA such that u £ C 1 ([0, 00), C 2JV+1 ). 
The dynamical system defined by 

S N (t) :u 0 £C 2N+1 -*u(t) £C 2N+1 , (59) 

possesses a bounded absorbing set Oq in <C 2N+1 and a global attractor An = w(0o) C Oq C C 2JV+1 . Moreover, 
for the absorbing ball, the entry time is independent of the initial data in C 2A,+1 , however large is the norm of 
the initial data. 

Proof: Let k = — m, m > 0. Taking the scalar product of m with u, we obtain the equation 

n=N n=N n=N n=N 

2 ffl ^2 K | 2 + A \{BdU)n\ 2 — 7 ^2 \u n \ 2 +m ^2 K| p+ 1 =0. (60) 

n=—N n=—N n=—N n=—N 


The operator — Ad : (C 2Ar+1 ,|| • H2) — > (C 2Ar+1 ,|| • 1 12) , defined by ©, is self-adjoint. If A* > 0 denotes the 
smallest eigenvalue for the eigenvalue problem —Adif = A*^, then since 


Aj = inf^ eC 2 iv+i 


(-A d 4’,ip)2 

(V’,'0)2 


= inf 


4>eC 2N + 1 


EL- 1 V IV’n+l - ifr, 


E n=N 
n=— 
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IV’n 




















































it is straightforward to check that if j < AJ, the dynamics are trivial in the sense limsup^^ ||it(i )||2 = 0. Thus 
we consider the case j > A*. From inequality GEJ we get that 

n=N 

ki\\u\\% +1 < m ^2 W P+ \ *i = mL i * £ , 

n=—N 


while from Young’s inequality we get that 

ki 


P~ 1 


7lMl a < T IMIS +po, n(N,i,p) = p + 1 Vtl(p + i) 


P 1 P+1 

7 p- 1 . 


Thus, m becomes 


1 d 


ki. 


2 dt M\i + ^\\n\r 2 + 1 <P0- 


(61) 


Now the result follows by applying Gronwall’s Lemma {22 Lemma 5.1, pg 167-168]: We get from Kill) that 


K*)Hi < 


2po 

ki 


( 2 


Now for any p\ satisfying 


V ki pt 


i 2 Po 

pl 1 'feT 


, for all t > 0 . 


(62) 


we derive from K>2I) . that for any set O of C 2Ar+1 , it holds S]y(t)0 C Oq for any 

/ , V _j_ \ -p 

2 / r, f 2p$ \ P+1 


t > to := — [Pi - -p 


hp 


ki 


where Oq is the ball of C 2W+1 of center 0 and radius pi. o. 

Remark 3.4 In the case of a DRGL equation with nonlinearity ( Q ) and k < 0, the functional E defined by (1431) 
is a Lyapunov function. Moreover, it can be checked (by following similar calculations as those of Proposition 
PHI) . that the set of stationary points £ is bounded. Hence, as it follows from 0E3I22], for each complete orbit 
containing uq lying in An , the limit set lo{uq) is a connected subset of £, on which E is constant. If £ is totally 
disconnected (in particular if £ is countable), any solution (Sjv(t)rto tends to an equilibrium point as t —► +oo. 
However, even in this finite dimensional problem for the DRGL, it could be an interesting issue, the investigation 
and rigorus justification of the bifuractions from the eigenvalue AJ, and convergence to (nontrivial) equilibria: 
writing the stationary DRGL problem as a real system, consisting of the (nonlinear algebraic) equations for the 
real part Re(u) and imaginary part Im(w), one could ask if the eigenvalue AJ could be a bifurcation point, from 
which two global branches bifurcate. These branches could consist of semitrivial solutions (i.e. solutions of the 
form (Re(n),0) or (0,lm(w)). Questions of this type will be considered elsewhere m- 


4 The case of an infinite lattice: Existence of global attractor for 
exponentially spatially localized solutions 

In this section, we prove the existence of a global attractor for the following complex lattice differential equation 

iu n + (& + ip){u„-\ - 2u n + Un+i) + (7 + iS)u n + if) + i()F(u n ) = g n , n £ Z, (63) 

Un(0) = tin, 0 , n G z. (64) 

We focus on the case of a gauge nonlinear interaction F(s) = |s| p- 1 s, p > 1. For specific values of the 
parameters, one recovers either the DCGL and DRGL equation or the weakly damped and driven DNLS. We 
refer to the pioneering work |S], on the existence of global attractors for lattice dynamical systems of first order 
and in 123 for extensions to various multidimensional lattices of first and second order. We remark that the 
discretization of the Laplacian, is not self-adjoint in a difference with [5] and the examples provided in E2, 
which increases considerably the manipulations needed, for the derivation of suitable estimates. The first result, 
is for the existence of an absorbing ball. 





Lemma 4.1 Assume condition (WS) on the weight function and that the parameters a, (3, 7 , 8, 77 , £ are chosen 
such that for some fixed e > 0 , 

a 0 := 6 - | - 20 - |d| DdT 1 ' 2 - \$\ f 1 + | > 0, C > 0. (65) 

Let (uo.n)nez = Uq £ Ig and (g n )ne% = g £ ig- A dynamical system can be defined by I6'dl) - f6'4l ), 

S(t) : uq £ il —> u(t) £ £$, ( 66 ) 

possessing a bounded absorbing set Bq in For every bounded set B of (f, there exists tn(B,Bo) such that for 
all t > t 0 {B, Bo), it holds S{t)B C Bo- 

Proof: We multiply © with O n u n , n £ Z add over all lattice sites, and keep the imaginary part. We obtain 
the equation 

— ^20n\u n \ 2 - ali(u„) - pl 2 (u n ) + 6^2 9 n \u n \ 2 + (^2 0 n \u n \ p+1 = 6 n g n u n , (67) 

nGZ nGZ nGZ nGZ 


where the terms Ii, I 2 , are defined as 

Il(^n) — ^ ^ ) n “h (-01/2 )n (-^^1^2 )n} — ^ ^ (^n+1 @n) (^l,n^2,n+l ^ , 2,n^’l,n+l) ? 

nGZ nGZ 

I 2 KO = E {(-Bui)„(.B6>m) w + ( Bu2)n(B0u2)n } = 2 E 6> ra (^ „ + u\ n ) 

nGZ 

~ EE"+! + d n ){u\ ,nUl,n+l + W2, n U2,n+l)- 


nGZ 


nGZ 


Using (WS), for the term U, we may get the estimate 

1-^1 (^n)| E ^ ^ |^n+l $n| |'^l,n'^2,n+l ^2,n^l,n+l| E D ^ ^ 6 n |'Wl,n'^2,n+l ^2,n^l,n+l | 
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nGZ 
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V.n(EZ 
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VnGZ 


r i/z i E 0 > 2 >™ 
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r 1/2 
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n |^2,n 

\n(=Z 

2 V , 2 +fx> 


\ 1/2 / \ 1/2 ' 

2 J ( E ^n|Ml,»+l| 2 j 

/ Vnez J 


1/2 
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^ ^ @n |^ 1 ,' 


^.nGZ 


Vn€EZ 


1/2 


< 2 M “ 1/2 E 0 ™l M bn | 2 E 0 


n|«2,n| 


VnGZ 

< AT 1 / 2 IMI? r 

For the second term on the rhs of I 2 , we have 


Vn£Z 


( 68 ) 


E(* n+l + 0n)(lll,nWl,ra+l + U2,nU2,n+l) 


nGZ 


< E $n+l |^l,nll'l,n+l | H” ^ ^ |'^l,n'^l,n+l | 


nGZ 


nG Z 


H"~ ^ ^ ^n+1 |'^2,n'^2,n+l | H"~ ^ ^ @n |^ i 2,n^ i 2,n+l 
nGZ n£Z 

/ \ !/2 / \ 1/2 

( E ^«+ii u i.™i 2 ) 


V.n(EZ 


E^+ll^l^+rl 2 ) + ^E/(El ft 
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1/2 


n |^1,n+l | 


+ I E! ^n+l| u 2,n 
\n£ Z 


i 2 ) 1/2 (e^+iI u2 .«+ii 2 ) 1/2 +(e 0 

/ VnGZ / \nGZ 


\ 1/2 / \ 1/2 

'n|^2,n| j [ E^l^ + l| j 

/ VneZ / 
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'n |^1,n+11 
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2 “I $n+l |lil, n |“ + El ^n+l|wi,n+l | 2 + El $n|wi jn | 2 + El 

LraGZ nGZ n£Z nGZ 

+ 21 J2 dn +^ U2 ^ 2 + X! ^n+l|'li2,n+l| 2 + E $n|^2,n| H - E ^n|^2,n+l| f 

LnGZ nGZ n£Z nGZ ) 

o ] ^E 0 n i Mi >"i 2 + E 6 ,?i + i i ui - ri + i i 2 + E 6 , "i Mi >"i 2 +^ _1 E ^ra+l|Ml,n+l | 2 


nGZ 


nGZ 


nGZ 


nGZ 


^ d ^ ^ |^2,n| + ^ ^ ^n+1 |^2,n+l | H - ^ ^ |^2,n| H” d. ^ ^ ^n+1 |'^ , 2,n+l | 

V nGZ nGZ nGZ nGZ 


, d d~ 

1+ 2 + ~ 

We insert (16811 . (1691) to 1671) . to get the inequality 

^IMI^ + 2c7o|M| 2 2 <-||.g|| 2 2. 

From GSJ> we derive that u G L°°([0, oo), £ g ): Gronwall’s Lemma implies that 

\\u(t)\\% < ||uo|||exp(- 2 a- 0 t) + -^—|| 5 |||{l-exp(- 2 f 7 0 t)}. 

6 6 Z(T q6 6 


(69) 


(70) 


(71) 


Letting t —+ oo we infer that 


limsup||t*(i)|| 2 2 < -- \\g\U. 

t —>oo 0 2 l<Toe 6 


Setting p 2 = ||g||^ 2 /2ero6 it follows that for any number p\> p the ball Bq of £ 2 g centered at 0 of radius p\ is an 

absorbing set for the semigroup Sit ): if B is a bounded set of £ g , included in a ball of £ g centered at 0 of radius 
R, then for t > to{B, Bo) where 

i R 2 


to = x— log ■ 


(72) 


2 cr 0 '~° p\- p 2 ’ 

it holds ||w(f )||? 2 < Pi, he. S(t)B C Bo- Note that in the absence of external excitation, the dynamical system 

exhibits trivial dynamics, in the sense that limsup t ^ ||u(f )||? 2 = 0 , as it follows from (17Tll . o 

£ e 

The next lemma provides us with the appropriate estimates, on the tail ends of solutions of m-m- 

Lemma 4.2 Let uq G B where B is a bounded set of £ g , and g G £ g . Moreover, we assume that the condition 
lt>5l on the parameters, is satisfied. Then, for any i) > 0, there exist T{if) and K(rj) such that the solution u of 
satisfies for all t > Tip), the estimate 


E Sn I Un W I 2 - f° r an y M > k (r)- 

\n\>2M a ° 

Proof: We consider a smooth function <f> G C' 1 (K + ,R), satisfying the following properties 

cf>(s) =0, 0 < s < 1 

0 < 4>(s) < 1, l<s<2 

<j>{s) = 1, s > 2. 


(73) 


and 


f> ( s ) |< Co, s G R + , 


(74) 


for some Cq G R. We shall use the shorthand notation <j) n = <j> We now multiply with the function 

<p n 9 n u n , nGZ, and we sum over all sites and keep the imaginary part. The resulting equation is 
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(75) 





where the terms Li, L 2 are found to be 
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Now for the second term of the rhs of we have 
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VnGZ / VnG Z / VnGZ / VnGZ / 


1/2 


< - E 0n+l^n+l |l^n | H - ^ ^ 0n+l^n+l l^n+11 H - ^2 4>nO n \u n \ 2 + X^+^n|'lt„+l| 2 > 

vnGZ nGZ nGZ nGZ J 

< Nd E 0n+l^n l^ n l "h ^ ^ 0n+l^n+l jl^n+l | H - ^ ^ 0n^n l^n [ H - d ^ ^ 0n^n+l |l^n+l | ^ 

v nGZ nGZ nGZ nGZ ) 



— ^ ^ l^n| H - ~ ^ ^ |^n| H“ ~ ^ ^ (0n+1 < / ) n)^n|'^?i 


n£Z 

4 _1 

2 


n£Z 


n£Z 


7-1 


H - ~ ^ ^ 0n+l^n+l |^n+l | ~ ^ ^ (0n+l 0n)^n+l |^n+l | 


< 


n£Z 

, d d _1 \ , „ , l2 c P i 

1 + ~ + -^- ) /2 &nOn\Un\ + “TT"- 


n£Z 


2 2 


n£Z 


M 


Inserting f7()ll and ITT!) to 117511 . we obtain the differential inequality 


1 d 

2 dt 


C , 1 


^ ^ f fin^n | | T<To ^ ^ (fn^n | U n | ^ M ^ ^ ^ 


n£Z n£Z 

Using Gronwall’s inequality, we obtain the estimate 


\n\>M 


(77) 


Y 4>n0‘ 


n un 


2 < e _2o ' o(t_4o) £ <M„ | Un (i 0 ) | 2 +^-[^pi + - E 0 " I 3 


n£Z 


n(zZ 


n | 2 | > 


|n|>M 


for t > to, where to is the time of entry of initial data bounded in £g, into the absorbing ball of radius p\ in Eg. 
Since g £ Eg, then for all 77 > 0, there exists K(i 7 ) such that 


9 /qi 1 

+ - E ^ I S" | 2 < V, VM > K{rj). 

\n\>M 


Therefore, for all g and for t > to and M > K(g), we obtain that 

E^ 0 " I Un I 2 - e~ 2a °^~ to) p\ + 2 -ry. 


n£Z 


Choosing t large enough, we may then obtain 


E 0™ | Un | 2 < E I Un I 2 - —' 

Z ' Z ' (Jf) 

\n\>2M neZ 


This estimate holds as long as t > T(g), where 


and M > K(g). This concludes the proof of the Lemma, o 

With Lemma IOI at hand, we are able to prove that the semigroup S(t) : (g —> Eg is asymptotically compact. 
The proof follows closely that of ;2J, adapted in the case of Eg, and is presented for the completeness of the 
presentation. 

Proposition 4.1 The semigroup S(t) is asymptotically compact in Eg, that is, if the sequence u n is bounded in 
Eg and t n —> 00 , then S(t n )u n is precompact in Eg. 

Proof: It follows from Lemma U~T1 that if u„ € Eg, such that || u n \\p g < r, r > 0, there exists T(r) > 0 and an 
integer N\(r), such that t n > T{r) for n > N\{r), and 

S(t n )u n C B 0 , Vn > N\{r). (78) 

From C3, there exists uo € £>o and a subsequence of S(t n )u n (not relabelled), such that S{t n )u n —*• uo in Bo. 
Lemma n~ 2 i implies the existence of some K\{g) and T\(rj) such that 

2 

E 0* I (S(t)S{T r )u n )i | 2 < y, t > Ti(r?), 

\i\>Kl(ri) 



where by (S(t)u n )i, we denote the *-th coordinate, of the infinite sequence S(t)u n G £g. Since t n —> oo, we may 
find N 2 (r, if) G N, such that t n >T r +T(rj), if n > N 2 {r, rj). Hence, 

E e t I (S{t n )u n )i | 2 = E e t I (S(t n - T r )s(T r )u n )i | 2 < 

|i|>^iW |i|>-K'i(*7) 

On the other hand, 

E * i m, r< 

\i\>K 2 (v) 

Choosing now K{rf) = m&x{K\{rf),K 2 {rj)), we get for all 77 > 0, that 

|| B{t n )u n Uq 11^2= ^ ( 0i | ( S(t n )u n U-o)i | T ^ | (S(tn)u n Uq)i | 

\i\<K(ri) \i\>K ( V ) 

<-y + 2 E 0 i(l {S{t n )u n )i | 2 + | (ti 0 )i | 2 ) < if 

\i\>K(r/) 

Note that the first estimate, comes from the strong convergence in the finite dimensional space C 2K ( r, ' +1 . o 
The main result of this section, which is a consequence of Proposition rm and El Theorem 1.1.1], can be 
stated as follows: 

Theorem 4.1 The semigroup S(t) associated to possesses a global attractor A = lo(B 0 ) C Bo C £g 

which is compact, connected and maximal among the functional invariant sets in £g . 

4.1 The finite dimensional approximation of the global attractor in £g 

This section, is devoted to the finite approximation of the global attractor, of exponentially localized solutions of 
m-m- Since from Theorems l2.2l and l4.ll the solution of m-m is in C 1 (K + , fg), an immediate consequence 
of the inclusion relation is that 

lim u n {t) = 0 , t > 0 . 

n —>oo 

Thus it is natural to seek for approximations of the global attractor, generated by the the following system of 
ordinary differential equations, supplemented with Dirichlet boundary conditions, 

iv n + (a + 0)(v n -1 - 2v n + i!„ + i) + (7 + iS)v n + {r) + i£)F{v n ) = g n , |n| < N, (79) 

w -(Af+l)(’) = «(2V+1) (0 = °> ( 80 ) 

u n (0) = u n , 0 , |n| < N. (81) 

System can be written as an evolution equation in C 2Ar+1 , this time endowed with the inner product and 
induced norm 

n=N n—N 

{u,v) 2e := Re E 9nU n v^, ||it|| 2e := E 0 n \u n \ 2 , u, v G C 2N+1 , (82) 

n=—N n=—N 

Since all the norms in the finite dimensional space C 2Ar+1 are equivalent, a result similar to Proposition rm 
can be produced, establishing the existence of global attractor in Xg := (C 2m+1 , || ■ || 2s ), with entry time 
independent of the initial data. However, since for the finite dimensional approximation, we are interested in 
a-priori bounds in C 2W+1 endowed with the || • || 2 „-norm, which should be independent of N, it is crucial to 
follow the procedure described in Lemma 1 m and pose the same conditions on the parameters. Thus working 
exactly as in Lemma rm we may prove the following 

Proposition 4.2 Let Vq := (tVi,o)|n|<iV G Xg. For 1 < p < 00 , there exists a unique solution of ^-M> 
such that v € C 1 ([0, 00 ), Xg). Assume further that condition < 65 1) holds. Then the dynamical system defined by 

S N (t) : ^0 G Xg —► v(t) G Xg, (83) 

possesses a bounded absorbing set Oq in Xg and a global attractor An = cj(Oq) C Oq C Xg : For every bounded 
set O of Xg, there exists t\{0, Oo) such that for all t > t\{0, Oo), it holds that SN{t)0 C Oo, and for every 
t > 0 Sn ft)A n = An- 


Following mumi (see also El EH for a similar idea applied to pdes considered in all of R^), we observe that the 
Absolution of can be extended naturally in the infinite dimensional space l 2 s , as 


{uN{t))N& 


v(t ) = (Un(t))|„|<JV, M < N, 

0, jnj > TV. 


(84) 


Let us note that in the light of JHS) the finite dimensional space Xg , could be viewed as a finite dimensional 
subspace of £#, with elements u £ £g satisfying the Dirichlet boundary conditions m - 

The global attractor A of the semigroup S(t) associated with (1651) - (1641 . will be approximated by the global 
attractor An of Sx{t) associated to 117^11 - (IKHIl . as N —> oo. Recall that the semidistance of two nonempty 
compact subsets of a metric space X, endowed with the metric dx(-, •), is defined as 


d(Bi,B 2 ) — sup inf dx(x,y). 

zeBi i/eft. 


Theorem 4.2 The global attractor An converges to A in the sense of the semidistance related to £g: we have 
that limAr^oo <1(An,A) = 0. 

Proof: We denote by U an open-neighborhood of the absorbing ball Bo of S(t). Obviously A attracts U. For 
arbitrary IV 6 N, we consider the semigroup Sx{t) defined by Proposition 14.21 and its global attractor An- 
Exactly as in Lemma f4. II it can be shown that Bo D Xg is also an absorbing set for Sjv(t). Therefore 

An C Bo n Xg c WD Xg, 

and An attracts U D Xg. In the light of Proposition 14.21 and [521 Theorem 11.2, pg. 28], it remains to verify 
that for every compact interval I of R + , 

5n{I) := sup supd(5jv(i)V ; 0) S(t)ipo) —>► 0, as N —> oo. (85) 

voGUDXe tei 

We consider the corresponding solution v(t) = S)v(f)vo, u(0) = vo, in Xg through (1791) - (IsHl) . Then by Proposition 
EH it follows that v(t) £ An for any t £ R + . Therefore, if p > 0 is the TV-independent radius of the absorbing 
ball On in Xg , then for every t £ R + , ||u(f)||| 9 < p 2 . Using (151) . we may construct the extension of v{t) in £g. 
The extension uxit) satisfies the estimates 

\\u N {t )\\ 2 e 2 < p 2 , ||wAt(t)||| < C(p,\\g\\q), (86) 

the latter derived by itTflll . According to ]52 Theorem 10.1 pg. 331-332] or [57] Lemma 4, pg. 60], for the 
justification of (1551) it suffices to show that UAr(f) converges to a solution u(t) of (1631) - (1641) in an arbitrary 
compact interval of R + , and uo = u{ 0) in a bounded set of £g. Let I be an arbitrary compact interval of R + . 
From estimates (1551) . we may extract a subsequence Uj of un, such that 

Uj{t) —^ u(t), in £g, as j —* oo , for every t £ D, (87) 

where D denotes a countable dense subset of I. 

For any t £ I we consider the sequence 

Xjv(t) := (u N (t),z) t 2, z £ £g, 

which by is differentiable as a function of t, and x'n( t) = («jv(t), %)e ■ Moreover, there exists £ € I such 
that, for fixed t, s £ I 

\Xn( t) — Xn (s)| = |(mjv (t) — un ( s ) ,z)( 2 \ = z)^\\t ~ s\ 

< supllujv^lUll^ll^lt-sl <C\t-s\. (88) 

?e/ e e 

i.e the sequence xn is equicontinuous. On the other hand it follows from that there exists an TV-independent 
constant C\ such that 


II u N (t) - Ujv(s )|| f 2 <Ci\t- s I 


(89) 


Hence by Ascoli’s Theorem, it follows that the convergence m, holds uniformly on I as N —> oo. Summarizing, 
we obtain for the subsequence Uj , the convergence relations 

Uj —> u in C(I, Cg), 

Uj ^u in L°°(I,^), (90) 

Uj —A u in L°°(I, ^). 

For the passage to the limit, we shall use an eqivalent formulation of m-m- Clearly, any solution of 
m - m , satisfies for every v £ £g and z(t) £ Cq°( I), the formula 


I ( iu(t ), v) ( 2z(t)dt + J ( Lu(t ), v) ( 2z(t)dt + J ((7 + i 8 )u(t), v) e 2z(t)dt 


+ J{(v + K)F(u{t)),v)qz(t)dt = J(g,v) e 2z(t)dt, (91) 

where F(u) = \u\ p ~ 1 u. Since for fixed N £ Z + , un is a solution of (t7§l) - (lfHTl . we may reproduce m for un, 
by multiplying (T 7 T 71 ) by v £ £g, in the fg-scalar product. By Lemma ITT! F : £g £g is Lipschitz continuous on 
bounded sets of l 2 g . Therefore from m, it follows that there exists a constant c(p), such that \ \F(uj)—F(u)\\(2 < 
c(p)\\uj — u\\t,2. Then from (19(71) . we infer 


(■ F (uj(t )) - F(u),v) t 2z{t)dt 


< 


\\F( Uj (t)) - F(u(t))\\(2 1 \v\\ t 2 Z {i)dt 


< C j || Uj{t) - u(t)\\ e 2 e \\v\\ e 2 g z(t) 


< cswp\\uj(t)-u(t)\\ e \\v\\ e 2 / \z(t)\dt-+ 0 , as /i —> 00. 
tei 6 e J 1 

On the other hand, we have from Lemma El that L : £g —> £g is globally Lipschitz on £ 2 e . Hence, we have that 


(Luj(t) — Lu(t),v)g2z(t)dt 


< 


\\Luj(t)-Lu(t)\\ e \\v\\<,2z(t)dt 


< L J Wujit) - u{t)\\ e 2 g \\v\\ t 2 g z{t) 

< Lsup \\uj(t) - u{t)\\ e 2\\v\\ e 2 / \z(t)\dt 0 , as /r- 

tei 9 6 J 1 


OO. 


Since / is arbitrary, iPill is satisfied for all t £ R + , i.e. u(t) solves (IfTTTl) - (IfT^l) . Moreover by (FTol) . we get that 
u(t) is bounded in £g for all t £ R+. Therefore u{t) £ A, which implies that Uj(0) — > m( 0), and it( 0 ) is at least, 
in a bounded set of £g. Since the convergence holds for any other subsequence having the above formulated 
properties, by a contradiction argument using uniqueness, we may deduce that the convergence holds for the 
original sequence un■ Condition m is proved, o. 

Remark 4.1 The global attractor of exponentially localized solutions for the DCGL and DNLS 
equations. To recover from the complex lattice differential equation (El, the DCGL equation 

Un - (A + ia)(u n -i - 2u n + u n+ i)+~iu n +[k + i(3)\u n \ p ~ 1 u n = fn, 1 < p < 00, 

M„(0) = Un,0, 


(with gauge nonlinear interaction (C/)-note that this time the nonlinearity is on the lhs of the equation ), we 
have to set d = a, $ = —A, 7 = 0, j = —7,7) = —/3, C = k and g n = if n . For these values of the parameters, 
condition El becomes 

^_l + 2 \-\a\Dd 1 / 2 -\(l + ^ + ^y k> 0 . 

In particular, in the case of the exponential weight 9 n = exp(At|n|) (exponential localization), we find d = e p , 
d = e -M , D = e p — 1. We set e = 2A > 0 and condition (1(551) reads as 


•7 > Xe p + 2|o| sinh(Ai/2), 7 < 0. 






Note that in the absence of the external excitation (/ n = 0), the dynamical system exhibits trivial dynamics, 
in the sense that limsup^g, ||w(t)||? 2 = 0. 

For the weakly damped and driven DNLS, (A = k = 0) 

u n -ia(u n _i - 2u n + u n+1 ) + iju n +(3\u n \ p ~ 1 u n = f n , nfZ, 1 < p < oo 

the condition on the dissipation parameter is 

—7 > 2\a\ sinh(/z/2), 7 < 0. 

Although this condition appears from the consideration of the problem in (f, it seems to be in conformity with 
the analysis of section 3.3, since there is not any restriction for the sign of the parameter (3. Let us note that 
even for the weakly damped and undriven NLS partial differential equation there exist initial data for which 
solutions may blow-up in finite time see EU This is never the case for the damped and undriven (/„ = 0) 
DNLS, for which solutions always exhibit energy decay. 
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